Is manganese-doped diamond a ferromagnetic semiconductor? 



Steven C. Erwin and C. Stephen Hellberg 
Center for Computational Materials Science, Naval Research Laboratory, Washington, D. C. 20375 

(Dated: February 2, 2008) 

We use density-functional theoretical methods to examine the recent prediction, based on a mean- 
field solution of the Zener model, that diamond doped by Mn (with spin S—5/2) would be a dilute 
magnetic semiconductor that remains ferromagnetic well above room temperature. Our findings 
suggest this to be unlikely, for four reasons: (1) substitutional Mn in diamond has a low-spin 
5=1/2 ground state; (2) the substitutional site is energetically unfavorable relative to the much 
larger "divacancy" site; 3) Mn in the divacancy site is an acceptor, but with only hyperdeep levels, 
and hence the holes are likely to remain localized; (4) the calculated Heisenberg couplings between 
Mn in nearby divacancy sites are two orders of magnitude smaller than for substitutional Mn in 
germanium. 

PACS numbers: 



I. INTRODUCTION 

A large class of dilute magnetic semiconductors (DMS) 
is based on manganese doping of III-V or group-IV hosts 
having the zincblende, wurzite, or diamond structure. 
Formally, Mn is either a single or double acceptor, de- 
pending on whether it substitutes for a group-Ill or -IV 
atom, respectively.^'^ Aside from this difference, there 
are no simple guidelines for predicting how the magnetic 
behavior of the resulting DMS depends on the choice of 
host semiconductor. Elucidating such guidelines would 
be of great practical interest for efforts to control and 
optimize the magnetic properties of DMS materials. In 
particular, it would be very helpful to understand how 
the resulting Curie temperatures depend on the choice of 
host semiconductor. 

An important and early contribution was made by 
Dietl, who used a mean-field solution of the Zener 
model to predict the Curie temperatures for Mn dop- 
ing of a wide range of II- VI, III-V, and group-IV host 
semiconductors!^ Mn ions with localized spin S=5/2 
were assumed to substitute on the cation site, and to in- 
teract by an indirect exchange interaction mediated by 
holes. Curie temperatures were calculated for fixed Mn 
content x = 5% and hole concentration p — 3.5 x 10^° 
cm~'^ using the mean-field result 

keTc = x^sNoSiS + 1)0^ Aprnkp /I2h'^ ■ (1) 

In this equation, the most direct dependence of Tc on 
the host semiconductor comes from the density per unit 
volume of cation sites iVg, which varies with the host lat- 
tice constant as A much weaker dependence on the 
host arises from the density-of-states effective mass m* 
and the Fermi wavevector kp. The remaining variables, 
i.e. the p-d exchange integral [3 and the Fermi-liquid pa- 
rameter Ap, were assumed not to depend on the host. 

The full numerical evaluation of Eq. (1) for different 
host semiconductors'* leads to values of Tc that indeed 
track the prefactor 1/a'^ quite closely.^ For example, the 
hosts with the largest lattice constants (CdTe, InAs, Ge) 
lead to the smallest predicted Curie temperatures (be- 



low '^75 K). Likewise, the hosts with the smallest lattice 
constants (ZnO, GaN, C) lead to the largest predicted 
values of Tc (well above room temperature). 

Diamond has the smallest lattice constant of all the 
semiconductors (a = 3.56 A) and, accordingly, the high- 
est predicted Curie temperature (Tc ~ 470 K). These 
extreme values make diamond a logical choice for exam- 
ining some of the assumptions underlying the predictions 
of high Curie temperatures for hosts having small lattice 
constants. In this paper, we use first-principles theoreti- 
cal methods to test the assumption that Mn occupies the 
substitutional site, with spin S = 5/2, in the diamond 
lattice. We find that substitutional Mn has a low-spin 
S=l/2 ground state in diamond. Moreover, because of 
diamond's small lattice constant, the substitutional site is 
energetically unfavorable relative to the much larger "di- 
vacancy site," denoted V2, which we identify as the more 
likely Mn impurity location. We analyze the electronic 
and magnetic structure of this divacancy Mn impurity in 
diamond, and show that it is unlikely to lead to magnetic 
ordering at any reasonable temperature. 




FIG. 1: (a) Ideal and (b) completely relaxed geometry of the 
Mn impurity in a diamond divacancy site. 



II. METHODS 

Most of the results reported in this article are for iso- 
lated Mn impurities, which we simulate using a supercell 
of 64 carbon atoms in the diamond structure. All atomic 
positions were relaxed using total energies and forces cal- 
culated within the generalized-gradient approximation 
to density- functional theory (DFT), as implemented in 
VASPiSii Carbon and Mn ultrasoft pseudopotentials were 
used with a cutoff energy of 286 eV. For total energies, 
we used 2x2x2 Monkhorst-Pack sampling of the Bril- 
louin zone, and 4x4x4 sampling for convergence checks. 
In the discussion of the electronic structure of the iso- 
lated Mn impurity at a divacancy site, eigenvalues at the 
zone center of a 128- atom supercell are reported in order 
that degeneracies are properly represented. 

Within the supercell formalism, the formation energy 
of a Mn impurity is given by 

[Mn"] ^ Et [Mn"] ~ nc - MMn + qEF; (2) 

where Et [Mn^] is the total energy of a supercell contain- 
ing nc carbon atoms and one Mn impurity, with chemical 
potentials /ic and ^Mn, q is the charge state of the Mn 
impurity, and Ep is the Fermi level. Since the host crys- 
tal is elemental, the chemical potentials are simply the 
energy per atom in the diamond phase of carbon and the 
ground-state aMn. 

The electrical activity of a Mn impurity is determined 
by its formation energy as a function of charge state. 
Electrically active defects will have more than one stable 
charge state within the host band gap; the value oi Ep 
for which two such charge states have equal formation 
energies is referred to as a "charge transition level." The 
charge transition level between the neutral and q — 
states is the donor ionization energy, while the transition 
level between neutral and q = —1 states is the accep- 
tor ionization energy. In practice, the total energy of a 
charged supercell must be calculated by adding a uniform 
compensating background charge to the supercell. The 
long-range nature of the Coulomb interaction then gives 
rise to a spurious Madelung-like contribution to the total 
energy, which must be subtracted from the calculated to- 
tal energy. This was done using a standard approach in 
which a multipole expansion of the defect charge (up to 
quadrupole order) was used to estimate the interaction 
energy analytically and then subtract it offi^ 

The magnetic interaction between Mn impurities in 
nearby divacancy sites was calculated using larger super- 
cells of 128 carbon atoms. No relaxations were performed 
for these calculations; tests for a few configurations con- 
firmed that this approximation did not change the results 
significantly. The interactions are represented by numeri- 
cally mapping the DFT total energies into the Heisenberg 
form. In practice this amounts to computing the differ- 
ence in total energy between the parallel and antiparallel 
spin alignment of two Mn impurities in a supercell. 
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FIG. 2; Formation energy of substitutional and divacancy Mn 
in diamond, as a function of constrained magnetic moment. 



III. ENERGETICS 

The Mn impurity enjoys a special status as a dopant 
in dilute magnetic semiconductors such as Mn^jGai-ajAs 
and Mn^jGei-a;. It preferentially occupies the substitu- 
tional (cation) site, where it serves two roles: it is electri- 
cally active as an acceptor, and it contributes a localized 
spin 5 = 5/2. At sufficient Mn concentrations, holes are 
created in the host valence band; these delocalized car- 
riers mediate the ferromagnetic interaction between the 
localized spinsi^ 

GaAs and Ge have similar lattice constants, and since 
strain contributes to the impurity formation energy, it is 
instructive to note that the local strain around a substi- 
tutional Mn impurity, denoted Musub, is very small for 
both hosts — the nearest neighbors distort by less than 
2% from their ideal positions. Tetrahedrally coordinated 
interstitial Mn impurities, denoted Mntot, have higher 
formation energies (by 1-2 eV for neutral impurities); 
interestingly, the atomic volume available at this site is 
the same as for the substitutional site, and consequently 
very little distortion is created by the Mn interstitial in 
either hostiS 

The lattice constant of diamond is 37% smaller than 
that of Ge, and so it is not surprising that a substitutional 
Mn impurity creates a much larger local strain in the 
lattice — within DFT the nearest-neighbor distortion is 
12%. Because of the extreme stiffness of diamond, such 
a large strain is expected to lead to a relatively high 
formation energy. A similar distortion occurs for Mn in 
the tctrahcdral interstitial site. 

A simple way to relieve this strain is to place the Mn 
impurity at a divacancy site — i.e. to remove two neigh- 
boring carbon atoms and place the Mn at their midpoint, 
as shown in Fig. 1(a). This doubles the atomic vol- 
ume available, leading to a much reduced local strain — 
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only 2% distortion of the nearest neighbors. Even in its 
ground-state configuration, which includes a slight off- 
center relaxation of the Mn ion as shown in Fig. 1(b), 
the local lattice distortion around this Mny^ impurity is 
very small. It is already widely believed that other large 
impurities — including Si, P, Ni, and Co — preferentially 
occupy the divacancy site in diamond jiStiiiiSiiiM hence 
our prediction that Mn also prefers this site appears plau- 
sible. 

To make the comparison between different sites quan- 
titative we turn now to the formation energies calculated 
within DFT. For neutral Mn impurities, we find the di- 
vacancy site is strongly preferred (by 0.6 eV) to the sub- 
stitutional site; the interstitial site is much less favorable 
than either (by ~10 eV). Figure 2 shows the formation 
energies for divacancy and substitutional Mn as a func- 
tion of magnetic moment, calculated using the fixed-spin- 
moment method within DFT. Even in the preferred di- 
vacancy site, the absolute formation is almost an order of 
magnitude larger for substitutional Mn in GaAs or Go, 
suggesting that standard methods for attaining high con- 
centrations of magnetic dopants will probably not work 
for diamond. 

The equilibrium magnetic moments AI in the substi- 
tutional and divacancy sites are also quite different com- 
pared to substitutional Mn in GaAs and Ge (M = 4 and 3 
IJ,B, respectively). In the substitutional site, Mn adopts 
a low-spin M = 1 state, with an extremely small en- 
ergy gain (0.1 eV) relative to the spin-unpolarized state. 



Hence, we anticipate that even if Mn could be forced into 
the substitutional site, its magnetic moment may not be 
sufhciently stable to facilitate magnetic ordering. In the 
divacancy site two nearly degenerate moments are found 
(A/ — 1 and 5 hb) with an energy difference of only 0.05 
eV; the formation energy vs. magnetic moment is sym- 
metrical about M — S ^B- This peculiar feature of the 
Mnva impurity will be explained in Section IV. 



IV. ELECTRONIC STRUCTURE 

The Mny^ impurity has an energy-level structure that 
arises from the relatively weak interaction of atomic Mn 
levels and host divacancy levels, as shown in Fig. 3. The 
host divacancy defect itself has been extensively stud- 
ied both theoretically and experimentally, and is well 
understoodiiS*iSiiLi^ The removal of two carbon atoms 
from the host lattice leads to a defect with D^d sym- 
metry, and thus the six available dangling-bond orbitals 
must form linear combinations that transform under the 
irreducible representations oig, a2u, e„, and For 
the neutral divacancy, six electrons are available to fill 
these levels. In an early theoretical study, Coulson and 
Larkins^ predicted the ground state one-electron config- 
uration to be a^gOj^e^. Later calculations confirmed this 
one-electron level ordering for both silicon and diamond, 
and we find it here too, as shown. 
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FIG. 4: Formation energy of divacancy Mn impurity in dia- 
mond, as a function of the Fermi level. Two charge-transition 
levels fall within the diamond band gap, as indicated. 



In Fig. 3 we consider separately the interaction of the 
majority and minority Mn levels with the four divacancy 
levels. Under the D^^ crystal field, the 3d level of Mn 
splits into one aig and two Cg levels. Hence, hybridiza- 
tion is allowed only between Mn and V2 levels both having 
either aig or eg symmetry, but no mixing is allowed be- 
tween Mn levels and V2 levels having either 0211 or sym- 
metry. Among the calculated majority Mny^ (t) levels we 
find very little mixing within any level, and therefore each 
has a clear parentage from either Mn or V2 orbitals, as 
indicated in Fig. 3. For the minority Mny^ (|) levels there 
is strong mixing between the empty V2 Cg level and the 
empty Mn eg(|) level, so that the occupied bonding com- 
bination falls well below the diamond valence-band max- 
imum (VBM). Considering for the moment only those 
levels that sit below the VBM, we find that seven elec- 
trons occupy majority levels and four electrons occupy 
minority levels. Hence, these low-lying levels carry a net 
"core" magnetic moment A/corc = 3/is. 

The highest occupied molecular orbital (HOMO) of the 
Mnv2 impurity has e„ symmetry and, because it can- 
not mix with Mn states, is purely derived from the V2 
dangling-bond orbitals. The two electrons in this four- 
fold level form a spin triplet (Mhomo = 1 Mb) shown, 
but the system is energetically indifferent to whether they 
are aligned parallel or antiparallel to the core spin (Fig. 3 
shows the parallel alignment; when the alignment is made 
antiparallel the spin-up and spin-down HOMO levels sim- 
ply switch positions). 

Because of this weak coupling between the HOMO 
spin and the core spin, the total magnetic moment of 
Mny^ can assume two possible low-energy values: either 

Mcorc + MuOMO = 5 /^s or Mcore " MrOMO = 1^-8 

for parallel or antiparallel alignment, respectively. More- 
over, there is an energy penalty for violating Hund's Rule 



by forcing the HOMO into a spin-singlet state; this is the 
origin of the energy increase that occurs when the total 
moment is constrained to 3 fj,B (the core moment alone). 
These two features — the weak core-HOMO spin coupling, 
and the Hund's Rule penalty for constraining the HOMO 
to a spin singlet — explain the unusual appearance, shown 
earlier in Fig. 2, of the Mny^ formation energy vs. total 
magnetic moment. 



V. ELECTRICAL ACTIVITY 

In Fig. 4 we show the calculated formation energies of 
Mny^ for charge states that are stable within the exper- 
imental 5.4-eV band gap of diamond. There are three 
such stable charge states {q = 0,-1, —2), making Mny^ 
formally a double acceptor. In principle this is similar to 
the case of substitutional Mn in Ge, which is also a dou- 
ble acceptor, and might suggest similar electrical activity 
for Mnv2 ■ 

Quantitatively, however, the two cases are quite dif- 
ferent. The first two acceptor ionization energies for Mn 
in Ge have been measured experimentally to be 160 and 
370 meV, relative to the valence-band edgeiSi For Mn 
concentrations of several percent, the width of the re- 
sulting impurity band is sufficient to create overlap with 
the valence band and thus to allow the holes to delocalize. 
Likewise, Mn in GaAs is a single acceptor with a mea- 
sured acceptor level at 113 meV,?^ which is accurately 
given by DFT calculations as 100 meV.^ 

On the other hand, the predicted first acceptor ioniza- 
tion energy for Mny^ in diamond is an order of magnitude 
larger, e(0/— ) = 1.2 eV, suggesting that the holes will re- 
main localized on the impurity sites for any reasonable 
Mn concentration. For this reason alone, it appears very 
unlikely that there will be sufficient carriers to mediate a 
ferromagnetic interaction between the localized spins of 
Mny^ impurities. 
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FIG. 5: Heisenberg spin-coupling constant between two Mn 
impurities in C (diamond). Si, and Ge host crystals. For 
diamond, Mn atoms are at divacancy sites; for Si and Ge, Mn 
atoms are on substitutional sites. 
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VI. MAGNETIC INTERACTIONS 

The weak coupling between core and HOMO spins in 
the isolated Mny^ impurity suggests that the effective 
coupling between spins of two impurities will also be 
weak. This expectation is confirmed in Fig. 5, where we 
show the results of our DFT calculations for (twice) the 
Heisenberg coupling constant between two Mny-j impu- 
rities in a 128-atom diamond supercell. For comparison, 
we show the results of similar calculations for substitu- 
tional Mn in Ge and Si. 

The results for substitutional Mn in Ge and Si show 
a similar trend. When the two Mn are nearest neigh- 
bors, there is a strong preference for antiparallel spin 
alignment. At larger separations the interactions vary in 
sign (Ge) and/or magnitude (Si), a consequence of the 
crystallographic anisotropy that arises from the strong d 
character of the acceptor wavefunction in substitutional 
Mn. The largest ferromagnetic coupling strengths are in 
the range 0.1—0.2 eV. (For Si there is a shift of all the 
couplings toward more favorable ferromagnetic interac- 
tions; since Mn-doped Ge is known to be ferromagnetic 
this suggests that Mn-doped Si might also be — perhaps 
with a higher Curie temperature). 

The coupling between two Mny-j impurities in dia- 
mond is substantially weaker: the largest ferromagnetic 
couplings are of order 1 meV, two orders of magnitude 
smaller than for Mn in Ge or Si. Hence, any ferromag- 
netically ordered phase involving Mny^ impurities would 
be characterized by a Curie temperature at most on this 
scale, i.e. of order 10 K. 
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VII. SUMMARY 

We have shown theoretically that the Mn impurity in 
diamond is energetically more favorable in the divacancy 
site than in the substitutional site. The magnetic prop- 
erties of the divacancy Mn impurity show nearly degen- 
erate low-spin and high-spin configurations; the degener- 
acy arises from the weak coupling of separate core and 
valence moments. Divacancy Mn is predicted to be a dou- 
ble acceptor with two hyperdeep levels, suggesting poor 
prospects for creating the delocalized holes in the valence 
band that are required to mediate a ferromagnetic inter- 
action within the Zener model. Explicit calculation of 
the effective Heisenberg coupling between Mn in nearby 
divacancy sites confirms that Mn-doped diamond shows 
little promise of ferromagnetism above a few Kelvins. 

Our findings may not be unique to diamond: other 
semiconductors with predicted high Curie temperatures 
have similarly small atomic volumes, which may likewise 
make the substitutional site unfavorable for Mn dop- 
ing. Comprehensive calculations of the site energetics 
for GaN, ZnO, and other host semiconductors would be 
of great interest. 
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